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Abstract.  A  direct  comparison  of  tlic  effectiveness  of  fission 
neutrons  at  high  d  1-0  3l-3rGy/min!  or  several  low  dose-rates 
(014  3-2cGy;min-  was  carried  out  under  identical  conditions. 
Monolayers  of  exponentially  growing  C3H/10T^  cells  were 
exposed  at  37  C  to  reactor-produced  neutrons  iflucnce-mran 
energy  £„  =  0-68MeY.  C  5"„  y  component,  frequency  mean 
linear  energy  _>F  =  21  keV, pm.  dose  mean  lineal  energy  )D 
=  42  ke\'//ini  in  an  8-pm  spherical  cavity  .  Survival  or  transfor¬ 
mation  induction  were  studied  at  five  doses  from  10-3  to  94  rGy. 
In  low  dose-rate  irradiations,  these  doses  were  protracted  over 
0-5.  1.  3  or  4-5  h.  resulting  in  17  dillcrcnt  dose-rates.  Up  to  six 
experiments  wrrr  performed  at  each  of  five  exposure  limes. 
Concurrently  with  transformation  we  studied  cell  proliferation  in 
control  versus  cells  irradiated  at  40cGy  iacutc  and  a  4-5-h 
protraction  and  found  no  evidence  of  a  shift  in  the  cell  cycle 
distribution  among  these  cells.  At  a  given  dose  and  dose-rate,  the 
cllcct  of  dose  protraction  on  survival  or  transformation  was 
assessed  by  the  dose-rate  modifying  (actor  ;1)RMF  .  defined  as 
the  low  :  high  dose-rate  rflect  ratio  at  the  same  dose.  Survival  or 
transformation  induction  curves  were  nearly  linear  with  initial 
slopes,  respectively,  of  about  fi-.ix  10'3  or  6-2  x  IO'6cGy''. 
Consistent  with  dose  response  curves.  I)R Mi  s  were  indepen¬ 
dent  of  the  dose  and  dose-rate.  The  mean  values  of  the  DRMF 
with  their  uncertainties  and  99"„  confidence  intervals,  based  on 
measurements  in  individual  doses  and  dose-rates  lor  survival  or 
transformation  were,  respectively:  101+003  0-92.  109  or 
0-98  +  0-0-1  0-83.  1-08  indicating  a  similar  precision  in  deter¬ 
mining  DR  Ml-'  for  survival  or  transformation,  and  no  dose  or 
dose-rate  influence  on  these  end  points. 


1.  Introduction 

Numerous  theoretical  and  experimental  studies 
have  been  conducted  with  the  objective  of  improving 
our  knowledge  of  dose  response  relationships  for 
radiation  carcinogenesis  at  high  dose-rate  or  pro¬ 
tracted  doses  of  medium  and  high  linear  energy 
transfer  (LET)  radiations.  Among  these  studies 
several  neutron  fields,  including  that  at  the  Armed 
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Forces  Radiobiology  Research  Institute  (AFRRI) 
TRIGA  reactor,  have  been  used  to  irradiate  C3H/ 
10Tj  ceils.  Exposure  of  these  cells  to  putative  carci¬ 
nogens,  including  ionizing  radiation,  may  convert  a 
cell  with  normal  growth  and  morphologic  character¬ 
istics  to  one  with  the  characteristics  of  a  tumour  cell,  a 
process  termed  neoplastic  transformation.  It  is  well 
established  experimentally  in  the  C3H/10j  system 
that  high-LET  radiations  have  greater  biological 
effectiveness  than  low  LET  radiations.  However,  the 
effects  of  protracted  or  split-dose  exposures  continue 
to  be  a  subject  of  controversy;  for  the  literature 
background  on  dose  protraction  studies  with  low  and 
high  LET  radiations,  see  Balccr-Kubiczek  and 
Harrison  1991c,  Brennerand  Hall  1990,  1992,  Elkind 
1991  and  Cao  el  al.  1992.  In  the  case  of  fission 
neutrons,  the  need  for  additional  transformation 
measurements  at  doses  <0-5 Gy  and  dose-rates 
$0-5cGv/min  has  been  identified  ;  Brenner  and  Hall 

1990,  1992,  Hall  1991,  Redpath  el  al.  1991). 

All  the  previous  dose-rate  experiments  with  C3H/ 
10TI  cells,  including  our  transformation  studies  with 
TRIGA  neutrons  (Balccr-Kubiczek  el  al.  1988,  1991, 
Balcer-Kubiczck  and  Harrison  1991a;,  used  a  con¬ 
stant  dose-rate  approach  in  which  the  total  dose  is 
proportional  to  the  exposure  time.  Consequently, 
different  irradiation  times  are  needed  to  accumulate 
different  doses,  a  procedure  which  makes  the  analysis 
of  the  protracted  data  conceptually  complicated 
especially  when  intervals  of  enhanced  sensitivity 
within  the  cell  cycle,  the  division  delay,  repair  or 
similar  kinetic  and  time-  or  dose-dependent  pheno¬ 
mena,  singly  or  in  combination,  arc  suspected  to 
affect  the  induction  process  (Rossi  and  Kellerer  1986, 
Brenner  and  Hall  1990.  1992,  Elkind  1991,  Hall 

1991.  Cao  el  al.  1992,  Watt  1992). 

The  objective  of  the  present  series  was  to  re¬ 
evaluate  the  effect  of  fission  neutron  dose-protraction 
on  neoplastic  transformation  of  C3H/10T^  cells  by 
using  a  constant  exposure  time  approach  in  which 
dose-rates  are  proportional  to  the  total  dose  (Kellerer 
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and  Rossi  1972).  Consequently,  different  doses  are 
accumulated  over  the  same  exposure  time.  An  addi¬ 
tional  advantage  of  the  constant  exposure  time  tech¬ 
nique  is  that  contributions  of  the  total  dose  and  dose- 
rate  can  be  separated,  so  consequently  questions  of 
time-,  dose-  and/or  dose-rate-dependent  windowed 
responses  can  be  addressed  directly  by  a  suitable 
choice  of  doses  and  exposure  durations  (Sykes  and 
Watts  1989).  In  our  present  experiments,  total  doses 
of  about  10,  25,  40,  60  and  90cGy  were  protracted 
over  t  =  0-5.  1,  3  or  4-5  h,  resulting  in  17  dose-rates 
from  014  to  3-2cGy/min.  Our  previous  constant 
exposure  time  experiments  with  X-rays  showed  that 
a  moderate  protraction  (r^5h)  of  doses  in  the  dose 
range  from  25  cGy  to  4  Gy,  can  produce  a  significant, 
statistically  discernable  reduction  of  cell  killing  and 
transformation  in  the  C3H/IOT4  assay  (Balcer- 
Kubiczek  et  al.  1987). 


2.  Materials  and  methods 

2.1.  Cell  culture  and  media 

C3H/10Ti  cells  from  the  line  established  by 
Rcznikolf  et  al.  1 1973a,  b)  were  maintained  as 
described  previously  ( Balcer-Kubiczek  et  al.  1988, 
Balcer-Kubiczek  and  Harrison  1990.  1991a  ci  in 
conformity  with  published  guidelines  for  this  assay 
(I  ARC/NCI /EPA  Working  Group  1985).  Our  qual¬ 
ity  control  of  experiments  with  C3H/I0TJ  cells  is 
described  in  the  reports  cited  above.  New  batches  or 
lots  of  scrum  arc  routinely  prescreened  using  pro¬ 
cedures  and  criteria  summarized  in  the  1  ARC/NCI  ’ 
EPA  Working  Group  (1985;  report.  Serum  lots  are 
characterized  as  being  suitable  for  the  transforma¬ 
tion  assay  by  ascertaining  compliance  with  survival 
and  transformation  expectations  after  4  Gy  X-rays, 
after  4  Gy  X-rays  with  post  exposure  to  01  /ig  TPA; 
ml.  or  after  2-5/rg  MCA/ml.  Results  for  two  lots.  - 
matched  with  respect  to  the  concentration  of  1 1 
natural  serum  components,  tested  for  use  in  the 
present  study  are  shown  in  Figure  1;  a  single  lot, 
denoted  scrum  B  in  Figure  1,  was  used  in  the 
experiments  reported  here. 


2.2.  Transformation  and  survival  assays 

Cells  in  passage  10  were  used  for  experimentation. 
Before  neutron  irradiation,  cells  were  plated  in  several 
25-cm2  flasks  at  5  x  I04  cells  per  flask.  Experiments 
were  performed  2  days  later  to  ensure  that  treated 
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Figure  1.  Comparison  of  transformation  induction  in  C3H/ 
10T}  cells  by  X-rays  or  ehrmirals  with  carcinogenic 
activity  (methylcholanthrcne.  MCA;  12-0-tetra- 
dcranoyl  phorbo!-!3-acetatc.  TPA’  to  demonstrate  the 
cflcct  of  serum  quality  and  concentration  on  transforma¬ 
tion  or  elonogenic  responses  idata  not  shown).  Other 
abbreviations  have  the  following  meaning:  10°„  +  5n„ 
FBS.  indicates  a  protocol  with  the  first  four  media 
changes  performed  using  medium  containing  I0P„  of 
foetal  bovine  scrum  iFBS  .  and  the  remaining  three 
media  changes  performed  using  a  5", ,-serum  rontent; 
10”..  FBS.  high-serum  medium-change  protocol. 


cells  were  in  exponential  growth  phase  (see  §  2-3) . 
Except  for  neutron  irradiation,  flasks  were  kept  in 
horizontal  position  at  all  times  to  facilitate  attach- 
ment/rcattachmcnt  of  mitotic  cells.  A  hypersensiti¬ 
vity  of  these  cells  to  transformation  was  suggested 
(Elkind  1991.  Brenner  and  Hall  1992,  C-ao  et  al. 
1992;.  Temperature  was  maintained  at  36  37°C,  in¬ 
cluding  during  neutron  irradiation.  Cells  were  trans¬ 
ported  to  the  reactor  in  a  20-litre  thermally  insulated 
container,  similar  to  that  described  previously 
(Balcer-Kubiczek  el  al.  1988,  Balcer-Kubiczek  and 
-  -Harrison  1991c!.  At  the  reactor  site,  cell  temperature 
before  and  after  neutron  irradiation  was  controlled 
by  immersing  culture  flasks  Idled  with  warm  medium 
in  a  constant-temperature  water  hath. 

Each  experimental  run  included  a  high  dose-rate 
exposurr  group  and,  in  most  cases,  more  than  one 
low  dose-rate  exposure  group  at  the  same  dose. 
Neutron  irradiations  were  arranged  such  that  low 
dose-rate  exposures  always  preceded  high  dose-rate 
exposure.  Consequently,  cultures  that  received  low 
dose-rate  irradiations  were  incubated  in  a  horizontal 
position  for  up  to  5h  before  transporting  them  from 
Bcthesda  to  our  cell  culture  laboratory  in  Baltimore. 
There,  cells  were  removed  from  flasks  by  trypsiniza- 
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Tabic  1.  Pooled  data  from  constant  exposure  time  experiments  with  TRIGA  fission  neutrons 


Mean 
dose  +  o 
(cGy) 

Mean  exposure 
time  ±  <7 
(min) 

Mean  dose 
ratc±<r 
(cGy/min) 

SF+o 

A’ 

Y 

X 

P 

T/.V+er 

(I0‘2) 

TR+o 

(io"*; 

10-68  +  0-74 

0-94  ±0-08 

1 1-35±0  77 

0-994±0-059 

692 

671 

23 

280 

3-310-7 

1-1010-24 

26-22+  1-32 

2-37  +  0-15 

1 1-0610-46 

0-780  +  0  1 76 

365 

342 

24 

297 

6-6±l-3 

2-19  +  0-46 

41-47+1-60 

1-39  ±0-08 

30-0212-34 

0-53810-063 

191 

165 

30 

334 

15-7±0-9 

4-1510-86 

62-28+1-84 

5-68  ±0-37 

11-0110-78 

0-43I±0-I00 

513 

453 

65 

259 

12-711-6 

5-6810-37 

9 1  -30  ±  4-45 

2-92  ±0-1 3 

31-291 1-28 

0-266±0-057 

333 

249 

89 

303 

26-712-8 

9-58  ±  1  -05 

10-23  ±  0-2 1 

29-87  +  0-19 

0-3410-0! 

0-94710-083 

293 

283 

10 

278 

3-411-1 

1-2510-39 

25-38  ±1-75 

29-34  ±1-09 

0-8710-06 

0-804  ±0-1 98 

272 

248 

24 

346 

8-811-8 

2-6710-55 

40-20  ±0-35 

28-16  +  2-30 

1-4410-11 

0-527±0-075 

359 

332 

27 

298 

7-5±l-5 

2-63  +  0-51 

61-43  +  2-17 

30- 10  ±0-07 

2-04±0-08 

0-546±0-067 

187 

159 

29 

245 

15-5±3-9 

6-6211-27 

91-95±  1-85 

28-84±l-16 

3-1910-06 

0-260  1  0-097 

76 

52 

24 

371 

31-616-5 

10-2412-10 

10  50  +  0-61 

57-56  +  2-18 

0  1810  01 

0-948±0-087 

311 

302 

10 

256 

3-2±l-0 

1  - 1510-38 

26-20+1-25 

57-83  +  2-18 

0-4510-01 

0-733±0-175 

228 

211 

18 

289 

7-911-9 

2-6810-65 

4l-50±  1-38 

6l-48±2-55 

0-68±0-04 

0-585  ±0-050 

264 

234 

33 

291 

12-512-2 

4-1510-76 

60-50  ±0-50 

59-1510-00 

1  -0210-01 

0-49910-119 

79 

67 

12 

333 

15-214-4 

4-95±l-43 

90-00  ±0-00 

62-53±2-15 

l-44±0-05 

0  341  ±0-051 

100 

67 

33 

360 

33-0±5-6 

11-1311-95 

25-92+  1-30 

179-3  ±6-6 

0-1410-0! 

0-806  ±0-1 45 

239 

224 

17 

307 

7-!±l-7 

2  1110-54 

40-53±0-9l 

186-4±4-8 

0-22±0-0! 

0-491±0-072 

221 

199 

31 

312 

14-0  ±2-5 

3-37±0-72 

61-34+  1-37 

178-3±!-0 

0-3410-01 

0-408±0-053 

125 

111 

17 

234 

13-613  3 

5-0811-36 

90-32  ±0-27 

180-6  ±4-3 

0-50±0-01 

0-232±0-093 

362 

291 

77 

246 

21-3±2-4 

8-8611-05 

40-40  ±0-60 

251-8±3-8 

0  1610-00 

0-661+0-079 

218 

203 

19 

263 

8-7  +  20 

2-71+0-70 

59-9.5  ±0-05 

258-8  ±  2-8 

0-23±0-0() 

0-432±0-008 

185 

158 

32 

280 

17-313-1 

5-6311-09 

90-45  ±0-90 

271-513-7 

0-3310-01 

0-250  ±0-022 

353 

279 

80 

299 

22-7±2-5 

7-8610-92 

Sl\  surviving  Iran  ion.  plating  ctlicicnry  of  control  was  55"„:  A'.  A”.  P.  total  number  of  dishes,  number  of  transformants  of  type  2 
or  3  'Reznikolf  tint.  1973b.  IARC, NCI/EPA  Working  Group  1985).  number  of  dishes  without  transformants,  number  of  rolls  per  dish 
corrected  for  plating  efficiency  of  unirradiated  cells,  or  surviving  fraction  and  plating  efficiency  of  irradiated  cells:  XI X,  transforma¬ 
tion  rate  per  dish  IARC  NCI/EPA  Working  Group  1985  and  its  standard  error,  o  (Hicbcr  el  at.  1987);  TK.  transformation  rate  per 
survivor  1  Han  and  Elkind  1979;  and  its  standard  rrror,  a  (Balccr-Kubiczck  el  at.  1987). 


lion,  counted,  diluted,  and  plated  at  concentrations 
estimated  to  result  in  either  250  viable  cells  per  dish 
for  the  transformation  assay,  or  50  viable  cells  per 
dish  for  the  cell  survival  assay.  As  in  previous  studies 
the  growth  medium  was  renewed  at  weekly  intervals 
(2  weeks  for  survival  assay  and  8  weeks  for  transfor¬ 
mation  assay,.  However,  in  the  present  experiments, 
medium  with  serum  content  reduced  from  10  to  5°„ 
was  used  in  the  last  three  medium  changes.  As  shown 
in  Figure  1,  this  modified  medium-change  protocol 
has  no  effect  on  transformation  frequency. 

Cell  survival  was  determined  by  colony  formation, 
while  neoplastically  transformed  foci  were  identified 
according  to  published  criteria  (Reznikoff  et  al. 
1973a,  b,  IARC/NCI/EPA  Working  Group  1985 i. 
The  end  point  of  transformants  per  surviving  cell  was 
calculated  by  the  null  method  of  Han  and  Elkind 
(1979)  with  uncertainties  determined  according  to 
our  modified  analysis  (Balcer-Kubiczek  el  al.  1987). 
In  addition,  we  determined  the  mean  number  of 
transformants  per  dish  as  recommended  by  IARC./ 


NCI/EPA  (1985)  with  uncertainties  calculated 
according  to  published  methods  (Hicbcr  et  al.  1987). 

2.3.  Proliferation  studies  of  irradiated  and  control  cells 

Concurrently  with  several  of  the  above  studies, 
experiments  were  performed  to  compare  prolifer¬ 
ation  of  cells  taken  to  and  from  the  reactor  facility 
(90  and  40cGy,  acute  or  a  4-5-h  protraction)  with 
those  continuously  maintained  in  our  laboratory. 
Up  to  four  cultures  per  condition  per  experiment, 
established  as  described  above  for  transformation 
assay,  were  used  to  establish  growth  curves  and  to 
perform  flow  cytometry  measurement  of  the  DNA 
content.  Growth  curves  were  initiated  by  plating  500 
cells  per  100-mm  dish.  Cell  counts  were  obtained 
every  12  h  for  a  period  of  2  weeks.  The  results, 
summarized  in  Table  2,  showed  no  difference  in 
growth  rates  and  in  the  cell  cycle  distributions 
between  travelled  and  non-travelled  cultures. 
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Talilr  2.  IVolifcTittion  parameters  of  C'.3H  lOTt  crib  mod  in 
reported  studies 


Paramrtrrt 

Travelled 

cells 

Laboratory 

control 

40cC,y 

brief 

40  flfiy 
-  4-5  li 

Cell  cvrlc  distribution  " 

G,i\l 

12 

10 

14 

12 

r.<i 

7-1 

70 

68 

s 

HI 

16 

16 

20 

Doubling 

lime  li 

IK  K  +  U  3 

Hi  2  ±U  2 

nd 

nH 

t  Cells  were  processed  immediatclv  after  return  from  the 
real  tor  fat  ilit\  ltd.  Not  done. 


2.1.  Ihisimrtn  and  irradiation 

Experiments  were  performed  at  tlic  Armed  Forces 
Radiobinlogy  Rescan  It  Institute  AFRRI  TRIGA 
Reactor  Facility  with  fission  neutron  field  optimized 
lor  low  dose-rate  studies.  The  technical  descriptions 
ol  the  neutron  source  and  ol’the  exposure  facility  were 
published  previously  Moore  and  Elsasser  198(5. 
Baleer-Kul)ie/ek  /■/«/'.  1988.  Red  pat  It  rial.  1991  .The 
mean  fluencr-weighlcd  neutron  energy  for  this  con¬ 
figuration  is  (H58  MeV  Goodman  1985  . 

Since  our  previous  experiments,  ilte  lead  shielding 
arrangement  around  the  point  ol  cellular  exposure 
was  redesigned  to  reduce  the  turnover  time  for 
samples  from  the  previous  If)  min  to  about  1  min.  and 
ret  I  nee  the  y  contamination  by  fully  enclosing  sam¬ 
ples  with  .Vent  lead  shielding  and  In  increasing  the 
souree-to-target  distance,  measured  from  the  central 
axis  of  the  reactor  core,  from  1  to  2  tn.  thus  providing 
favourable  conditions  lor  the  multiple  protracted 
exposures  required  in  the  present  study.  The  latter 
modification,  suggested  by  one  of  us  (i.  H.H.  . 
permits  achieving  dose-rates  of  <()•*> c(»y  min  at 
correspondingly  higher  reactor  power,  leading  to  an 
improved  neutron-io-total  dose  ratio. 

Dosimetric  measurements  were  performed  at  the 
plane  of  cellular  exposure  prior  to  each  of  our  Iff 
experimental  runs  using  the  paired  ionization 
chamber  technique  ICRl  197b,  AAFM  1980, 
Goodman  198a.  /.email  rl  al.  1988  .  A  lissuc- 
eciuivalent  I  K  -wall,  methane-based  TK-gas  filled 
0-5-cnr'  Kxradin  chamber  was  used  to  measure  y  and 
neutron  events,  whereas  a  0-5-cm3  magnesium-wall, 
argon-filled  chamber  was  used  to  measure  the  y 
events  in  the  neutron-plus-y  field.  From  the  measure¬ 
ments  made  by  the  paired  chambers,  the  portions  of 
neutrons  and  y  in  the  mixed  field  could  be  deter¬ 
mined.  For  the  present  irradiations,  the  ratio  of 
neutron  to  total  dose  was  0-95+  1  at  all  dose- 
rates.  The  instrumental  uncertainties  for  dose  and 


dose-rate  determinations  were  +3-1  and  3-G"„. 
respectively. 

Recently,  the  radiation  quality  of  the  TRIGA- 
rcactor  neutron  field  has  been  characterized  in  terms 
of  microdosimetric  spectra  (H.  Gerstenberg.  manu¬ 
script  in  preparation'.  The  spectrum  was  calculated 
from  a  TRIGA  neutron  spectrum  based  on  both 
measurements  and  calculations  Aerbinski  et  al. 
1981a.  hi.  Briefly,  the  mean  values  of  lineal  energy 
based  on  measured  spectra  using  I-  or  8-/<m  dia¬ 
meter  detectors,  respectively,  are:  ty  =  38-9  kcY//tm 
and  vD  =  70-5  kcV/yim,  or  vf  =  21l  kr\’//mi  and 
vD  =  41-5keY//tm. 


2.5.  Statistical  analysis 

At  the  total  doses  used  in  this  study  the  survival 
level  were  >  20°, ,.  Accordingly,  a  linear  quadratic 
equation  —  ln.V=at  ■d+fi  d"‘  was  appropriate  to 
analyse  survival  data.  In  the  case  of  transformation 
data,  several  radiobiologically  realistic  relationships 
were  tried,  such  as: 


T=a  +  b  -  d 

(1) 

T=  a  +  b  ■  d  + 1  ■  d2 

(2) 

T=hd+cd2. 

13) 

where  T  is  the  incidence  of  transformation,  d  the 
dose.  a.  b  and  r  are  regression  coefficients  for  the 
induction  of  neoplastic  transformation.  The  calcu¬ 
lated  best-fit  curves  were  compared  using  several 
goodness-of-fit  criteria.  The  coefficient  of  determina¬ 
tion  was  used  to  estimate  the  fraction  of  the  total 
variance  accounted  by  a  model,  and  the  model 
selection  criterion  Akaike  I97»>  was  used  to  evalu¬ 
ate  least -squares  fittings  to  competing  models.  Cal¬ 
culations  were  performed  using  commercial  software 
PSI-Plot.  Poly  Software  Inc..  1993  . 

Quantitativ  e  comparison  of  dose  response  curves 
liir  cell  transformation  from  our  five  protraction 
protocols  was  performed  by  multiple-regression 
analysis.  We  used  this  procedure  to  lest  specifically 
whether  the  high  dose-rate  regression  line  and  pro¬ 
tracted  irradiation  regression  lines  could  have  come 
from  populations  with  the  same  regression  coeffi¬ 
cients  Klcinbatim  et  al.  1988.  Peixoto  1993).  Com¬ 
mercial  software  MathCAD  version  2-5,  MatliSoft, 
I uc.,  1989  was  used  for  interactive  programming  of 
/■-value  calculations  from  transformation  data  sets  in 
Table  I  'see  Appendix  . 

In  addition,  the  effect  of  dose  protraction  on 
survival  or  transformation  was  assessed  by  the  dose- 
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NEUTRON  DOSE (cGy) 

Figure  2.  Survival  of  C3H  1011  rolls  after  irradiation  with 
AFRR1  TRKJA  neutrons.  Symbols  correspond  to  the 
following  patterns  of  exposure:  Q.  acute:  and  •. 

t  =  0-51i:  A.  t  =  I  h:  T.  r  =  3h:  and  ♦  .  r  =  4-5 h:  t 
denotes  irradiation  times  with  the  actual  values  listed  in 
Table  1 


rate  modifying  factors  'DRMF).  defined  as  the 
low :  high  dose-rate  effect  ratio  at  the  same  dose. 
Standard  formulas  were  used  to  calculate  the  uncer¬ 
tainties  in  DRMF  estimates  shown  in  Figures  4  and 
5.  The  DRMF  for  survival  or  transformation  were 
analysed  by  linear  regression  to  determine  whether 
or  not  they  were  independent  of  dose  and  dose-rate. 
The  actual  equations  were: 

DRMF—  a  +  h-d  4 

I)RMF=a  +  b  -5  ;5 

where  d  the  dose.  S  the  dose  rate,  a,  b  are  regression 
coefficients  in  lable  3.  The  /-values  obtained  from 
data  analysis  by  the  above  equations  1  and  5 
were  compared  with  the  tabulated  values,  as 
described  in  the  Appendix.  Required  calculations 
were  performed  using  commercial  software  Math- 
CAD  version  2-5,  MathSoft.  Inc..  1989  . 

3.  Results 

As  shown  in  Fable  1  and  Figure  2.  there  were  no 
significant  effects  of  temporal  dose  distributions 
on  cell  survival.  The  shape  of  the  curve  was 


close  to  exponential  with  a  linear  coefficient  a 
=  (6-48+  1-46)  x  10-3cGy_I  and  a  quadratic  coef¬ 
ficient  (1-03  +  0-25)  x  I0~4cGy”2.  These  results 
are  in  good  agreement  with  observations  from 
related  experiments  with  high-LET  radiations 
reported  earlier  by  us  (Balcer-Kubiczek  et  al.  1991, 
1993,  Balcer-Kubiczek  and  Harrison  1991a)  as  well 
as  others  (Han  and  Elkind  1979,  Hill  et  al.  1982, 
Hieber  et  al.  1987,  Saran  el  al.  1991,  Goodhead  et  al. 
1992,  Komatsu  et  al.  1993). 

As  in  the  earlier  studies  with  cycling  or  stationary 
cultures  of  C3H/10T^  cells  and  TRIGA  neutrons 
(Balcer-Kubiczek  et  al.  1988,  and  references  cited 
above),  the  dose-response  curve  for  transformation 
induction  was  well  represented  by  two-  or  three- 
parameter  models,  given  by  the  equations  (1 ),  (2)  or 
(3),  whether  exposure-time  groups  were  considered 
separately  or  pooled.  However,  the  amount  of  data 
in  each  dose  response  set  was  too  small  to  assert  any 
model  conclusively  by  the  usual  curve  fitting. 
Multiple-regression  analysis  of  all  the  transforma¬ 
tion  data  in  Table  1  showed  F— 1-59,  whereas 
^o-95  16)  — 3*01  and  F0.gq'A.  1 6  i  =  5-24  are 


'si- 


Figuiv  3.  Neoplastic  transformation  of  C3H/I0T)  cells  after 
irradiation  with  AFRRI  TRKJA  neutrons.  Symbols 
correspond  to  the  following  patterns  of  exposure:  O, 
acute:  •,  t  =  0-5:  A.  T=lh:  T,  r  =  3h;  and  ♦, 
T  =  4-3h:  T  denotes  irradiation  times  with  the  actual 
values  listed  in  Table  1. 


BEST  AVAILABLE  COPY 


561 


E.  K.  Baker -Kubicztk  et  at. 


Figure  4.  Dose-rale  modifying  factors  for  survival  of 
( "tl  I  IDTJ  crib  plotted  as  a  fuiirlion  of  AFRRl  'I'RlfiA 
neutrons  d< >se  rate.  Symbols  correspond  |o  the  following 
patterns  of  exposure:  •.  T  =  0  .1  li:  A.  t  =  I  h:  ▼ .  1  =  3  lr. 
and  ♦.  i  =  A -.Mi:  denotes  irradiation  times  with  the 
actual  values  listed  in  Table  I. 


required  lot  differences  among  the  means  to  be 
significant  at  the  5  and  1 levels,  respectively.  Thus, 
we  concluded  that  the  live  dose  response  curves 
were  satisf'aetorilx  homogeneous,  and.  consequently, 
that  at  a  given  dose  there  was  no  dependence  of 
transformation  response  on  dose  protraction. 

Accordingly,  we  were  able  to  pool  the  data  from 
out  live  exposure  time  groups  in  determining  dose 
dependencies  in  transformation  induction.  For  the 
models  1  3  .  the  goodness-of-fit  measures  had 

similar  values:  coefficients  of  determination  003 
0  <♦:».  model  selection  criteria  2-5  2-8.  correlation 
coellie irnts  001  0  05.  I’rcdietahlv.  the  numerical 
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figure  ").  Dose-rale  tnndilx ing  factors  lot  the  induction  of 
neoplastic  tr.iMsIonn.itiiMi  til  ("111  |0  |"l  eell'  plotted  as  .i 
ftiueiioii  of  AFRRl  TRKiA  neutrons  dose  rate.  Svmbob 
eorrespoud  to  the  following  patterns  of  exposure:  #. 
T  =  11  alt:  A.  t  =  I  It:  ▼.  r  =  3h:  and  ♦.  r  =  4:>h: 

denotes  irradiation  times  with  the  actual  values  listed  in 
Table  I. 


Table  3.  Linear  repression  eoellieients  and  the  test  of  indepen¬ 
dence  for  dose  rate  modifying  factors  DR  Ml 


a 

bt 

/ 

Survival 

DRMF  versus 

dose 

(1-95.5 

IHM1I 

023 

DRMF  versus 

dose-rate 

0  975 

0  071 

(Mill 

Transformation 

DRMF  versus 

dose 

(i'll.! 

-7  x  111  4 

(HR 

DRMF  versus 

dose-rate 

HH4 

Otl.V.' 

or. 

a.  intrreept:  b.  slope:  /  value  delined  in  Appendix, 
t  The  unit  to  be  used  with  factor  b  is:  lor  the  dose  response 
relationships  e(i\  “  1  and  for  the  dose-rate  response  relation¬ 
ships  :min'eOy  .  Dose-rate  modifving  laetors  calculated  from 
the  data  on  survival  or  transformation  Table  I  are  plotted, 
respectively,  in  Figures  4  nr  5  as  a  him  lion  of  dose-rate. 


values  of  regression  coefficients  depended  on  the 
model.  Our  best  estimate  of  the  initial  slope  based  on 
the  ponied  data  is  50  7-5x  lO'^rGy"1,  and  the 
mean  weighted  by  the  reciprocal  of  the  variance  is 
6-2  x  10’^cGy"1.  The  solid  line  in  Figure  3  repre¬ 
sents  the  best  fit  to  the  linear  quadratic  equation  (31, 
with  the  numerical  solution  b  =  7-5  +  2-li  x  10'h 
cGy-\  and  r=  f2-4±  1-9  10'HcGy‘2.  Thus,  if  the 
statistical  analysis  is  restricted  to  the  initial  part  of 
transformation  induction  nine  <60cGy>  a  linear 
fit  is  an  adequate  representation,  since  the  contribu¬ 
tion  from  the  quadratic  term  is  <  10",,. 

To  further  confirm  the  above  findings,  dose-rale 
modifying  factors  for  survival  or  transformation, 
shown  in  Figures  4  or  5.  were  tested  lor  dependence 
on  dose  and  dose-rate  by  the  equations  { 4 '  and  (5). 
Specifically,  the  dependence  of  DRNIF  on  dose 
provides  qualitative  information  about  the  basic 
form  of  dose  response  curve  obtained  for  different 
dose  protraction  patterns,  relative  to  the  basic  Ibrm 
of  (lose  response  curve  obtained  from  a  high  dose- 
rale  exposure  at  the  same  dose.  The  results  of  the 
numerical  curve  lining  of  our  lour  data  sets  to  a 
straight  line  with  the  intercept  a  and  the  slope  h 
cGy~  1  or  min/cGy  are  shown  in  Table  3.  The  last 
column  in  'Fable  3  lists  /-values  of  about  0-2  for 
dependence  of  DR  MFs  on  the  dose  and  ofahottt  0-5 
for  dependence  of  I)RMFs  on  dose-rate.  At  3  =  0  05, 
the  eritiral  /, 


1)015 


15'  =  2131  and  at  3  =  001, 
/<)(»(> <  15  =2-947.  Therefore,  we  are  required  to 
conclude  that,  at  the  5  or  i"„  levels  of  .significance, 
the  1)R  MFs  were  independent  of  both  the  dose-rate 
and  dose:  that  is,  the  basic  forms  of  dose  response 
curves  are  the  same.  Intuitively,  the  same  conclusion 
can  be  reached  by  observing  that  both  the  intercepts 
and  the  mean  values  of  DR  MFs  are  close  to  unity; 
for  survival  and  transformation  ihev  were  respect¬ 
ively:  101  (0  02.  109:  and  0  08  (0  83.  I -08),  where 
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the  numbers  in  brackets  are  99%  confidence  limits. 
By  inspecting  Figures  4  and  5,  an  additional  obser¬ 
vation  can  be  made  that  the  quality  of  the  transfor¬ 
mation  data  can  be  quite  good,  comparable  with  the 
quality  of  survival  data.  Thus,  the  scatter  in  transfor¬ 
mation  data  sets  should  not  be  relied  upon  in 
justifying  a  wide  range  of  model-based  speculations 
(Brenner  and  Hall  1990,  1992,  Hall  1991). 


4.  Discussion 

In  this  report  cell  survival  and  the  induction  of 
neoplastic  transformation  were  studied  after  fission 
neutron  irradiations  at  up  to  five  exposure  durations 
per  dose.  Such  an  experimental  design,  the  first 
reported  for  high-LET  radiation  studies,  permitted 
direct  examination  of  the  dose-rate  influence  on 
these  end  points  at  multiple  dose  rates.  The  dose-rate 
range  included  single  values  of  dose  rate 
i$0-5cGy/min  previously  studied  in  this  and  other 
laboratories  using  the  constant  dose-rate  approach. 
The  reference  acute  irradiations  were  performed 
concurrently  with  protracted  irradiations  at  the 
same  dose.  Identical  protocols  were  used  in  all 
irradiations,  including  media  composition,  age  and 
proliferative  state  of  cultures,  initial  plating  densities 
and  the  numbers  of  dishes,  thus  providing  statisti¬ 
cally  acceptable  conditions  for  testing  the  null 
hypothesis  that,  at  the  same  dose,  there  was  no 
difference  in  the  effectiveness  of  fission  neutrons  due 
to  differences  in  the  temporal  pattern  of  exposure. 

Technical  limitations  in  the  present  study  pre¬ 
cluded  investigations  of  neutron  doses  of  <  lOcGy,  a 
limitation  that  also  applies  to  most  (Hill  et  al.  1982. 
1984.  Balcer-Kubic/ek  el  al.  1988,  Goggle  1988, 
Miller  et  al.  1990,  1993.  Hill  and  Zhu  1991,  Miller 
and  Hall  1991.  Sarane/a/.  1991)  but  not  all  in  vivo  and 
in  vitro  studies  of  radiation  carcinogenesis  (Ullrich 
1984,  Balcer-Kubiczek  and  Harrison  1991a).  Never¬ 
theless,  two  interesting  and  potentially  important 
results  emerged  from  the  present  investigations.  First, 
we  confirmed  that  the  dose-response  curves  for 
survival  and  transformation  at  doses  <  1  Gy  are  best 
described  by  a  linear  equation  or,  possibly,  a  linear 
quadratic  equation  with  a  weak  positive  quadratic 
term.  Similar  to  previous  results  with  exponential  or 
plateau  cultures  of  C3H/10T^  cells  and  TRIGA 
neutrons  (Balcer-Kubiczek  and  Harrison  1991a)  as 
well  as  with  some  other  neutron  sources  (Coggle 
1988,  Saran  et  al.  1991,  Komatsu  et  al.  1993),  we  did 
not  observe  ‘supralinearity'  in  the  dose-response 
curve  for  transformation,  often  associated  with  the  so- 
called  ‘inverse  dose-rate  effect’  reported  Tory-,  X-rays 
(Borek  and  Hal!  1974,  Miller  and  Hall  1978.  Little 


1979,  Hall  and  Miller  1981)  and  fission  neutrons 
(Hill  el  al.  1982,  Ullrich  1984).  Indeed,  models  by 
Rossi  (1981),  Rossi  and  Kellerer  (1986),  Burch  and 
Chesters  (1986),  Brenner  and  Hall  (1990)  and  Elkind 
(1991)  show,  either  directly  or  by  implication,  that 
the  enhancement  of  an  effect  due  to  protraction  may 
occur  in  a  dose  range  w'here  a  dose-response  curve 
has  a  fine  structure  such,  for  example,  has  a  locally 
negative  (convex)  curvature;  ‘bending  down’  seg¬ 
ments  of  response  curves  after  brief  and/or  protracted 
exposures  to  low  and  high  LET  radiations  were 
reported  for  neoplastic  transformation,  animal  carci¬ 
nogenesis  and,  more  recently,  for  mutation  induc¬ 
tion,  but  not  cell  survival  or  mortality  (Miller  and 
Hall  1978,  1991,  Little  1979,  Hall  and  Miller  1981, 
Hill  et  al.  1982,  Ullrich  1984,  Bettega  el  al.  1992, 
Tauchi  et  al.  1 993 ) .  Thus,  our  second  important  result 
that,  at  a  given  dose,  survival  or  transformation  did 
not  vary  with  a  temporal  pattern  of  exposure,  can  be 
seen  as  a  direct  consequence  of  the  first  result,  just 
discussed  above,  that  dose-response  curves  for  neu¬ 
trons  are  linear  <  1  Gy.  Similar  conclusions  were 
published  by  Coggle  (1988)  based  on  the  analysis  of 
rodent  carcinogenesis  data. 

Recently,  several  interpretations  of  locally  com¬ 
plex  dose -response  curves  and  the  inverse  dose-rate 
effect  have  been  proposed  (Rossi  and  Kellerer  1986, 
Dennis  and  Dennis  1988,  Sykes  and  Watt  1989, 
Brenner  and  Hall  1990,  Elkind  1991,  Watt  1992). 
Most  of  these  models  are  an  extension  of  an  idea 
(Oftedal  1968)  that  there  is  a  time  period  in  the  cell 
cycle  during  which  cells  are  more  susceptible  to 
transformation  (or  mutation;  Tauchi  et  al.  1993)  than 
al  other  times.  A  corollary  is  that  for  exponentially 
growing  cultures  with  the  same  distribution  of  cells 
throughout  the  cell  cycle  at  the  beginning  of  expo¬ 
sure,  different  distributions  and,  therefore,  different 
fractions  of  the  hypersensitive  cells  arc  expected  at 
the  end  of  protracted  versus  brief  exposures.  At 
present,  there  is  no  agreement  among  the  models  as  to 
a  duration  (Rossi  and  Kellerer  1986,  Brenner  and 
Hal!  1990,  Elkind  1991,  Harrison  and  Balcer- 
Kubiczek  1992,  Watt  1992)  or  a  specific  location  of 
the  sensitive  period  within  the  cel!  cycle  (Brenner  and 
Hall  1990,  Elkind  1991,  Cao  el  al.  1992,  Komatsu 
el  al.  1993,  Miller  et  al.  1993,  Tauchi  el  al.  1993),  but 
periods  of  sensitivity  lasting  from  a  few  seconds  or  less 
to  about  1  or  2  h  can  be  inferred,  as  discussed  below. 
Although  the  present  study  was  concerned  with 
investigation  of  dose-rate  effects  of  fission  neutrons 
rather  than  w-ith  testing  the  biophysical  models  cited 
above,  our  measurements  of  the  cell  cycle  distribu¬ 
tions  of  control  and  irradiated  cells  (Table  2)  argue 
directly  against  the  validity  of  these  notions.  Inde¬ 
pendently,  quantitatively  similar  results  and  conclu- 
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sions.  based  on  experiments  with  fission  neutrons 
(Saran  et  at.  1 993 1  and  with  monocnergctic  neutrons 
of  several  different  energies  (O  f).  1-0  and  6-OMcYi 
(Pazzaglia  cl  al.  199.3,.  were  recently  reported  for 
control  and  neutron-irradiated  C3H/I0T]  cells. 

In  experimental  investigations  of  dose-rate  or 
dose-fractionation  effects,  it  is  important  to  ascertain 
whether  the  variability  of  a  response  ran  be  attri¬ 
buted  to  a  pattern  of  exposure  and  not  merely  the 
scatter  of  data.  Figures  -I  and  5  show  enhancements 
(and  reductions  of  less  than  a  factor  of  2  at  several 
dose-rates,  such  as  are  reported  by  some  investiga¬ 
tors  as  the  inverse  dose-rate  effect.  However,  in 
contrast  with  our  present  report,  these  other  con¬ 
stant  dose-rate  studies  were  limited  to.  at  best,  two 
dose-rates  'Hill  dal.  1982.  Hill  and  Zhtt  1991.  Miller 
and  Hall  1991.  Redpath  et  al.  1991  and  or  a  narrow 
range  of  doses,  including  measurements  at  a  single 
dose  I’llrieh  198  !.  Hill  and  Zhu  1991.  Miller  rial. 
1991.  1993.  Redpath  rial.  1991  :  consequently,  these 
results  do  not  lend  themselves  to  rigorous  analysis  by 
conventional  statistical  methods  such  as  presented 
here,  nor  to  the  analysis  by  multiparameter  models, 
such  as  those  proposed  for  transformation  data  by 
Burch  and  Chesters  1986  .  Rossi  and  Kellercr 
1986  or  Brenner  and  Hall  1990  . 

We  examined  the  latter  model  analytically  Harri¬ 
son  and  Baleer-Kubiezek  1992  with  the  intention  of 
applving  the  model  to  our  transformation  data  lor 
IRK  i.\  neutrons  Baleer-Kubiczck  elal.  1988.  1991. 
Balcer-Kubiczek  and  Harrison  1991a.  and  thi- 
report  and  found  a  lack  of  a  single-valued  solution 
fin  the  duration  of  the  postulated  transformation- 
susceptible  portion  ol  the  cell  cycle  lot  neutrons 
Figures  I  and  2:  Harrison  and  Balcer-Kubiczek 
1992  :  t hi--  is  in  contrast  with  conclusions  published 
by  these  authors  that  the  model  predicts  a  single 
sensitivity  period  of  I  -02  +  0-25  h.  Our  subsequent 
examination  of  the  mathematical  determination  of 
model  input  parameters  from  experimental  data 
Balcer-Kubiczek  and  Harrison  1992  rescaled  an 
additional  caveat,  that  the  prediction  of  the  inverse 
dose-rate  for  transformation  by  neutrons  requires  ;m 
unrealistically  low  value  of  frequenev  mean  lineal 
energy,  such  as.  q  =  1  keV  /on  in  an  8-//m  diametet 
cavity  fin  fission  neutrons.  As  shown  here,  this 
assumption  is  out  of  line  for  TRKiA  neutrons  with 
_q=2l-2kc\  pm.  In  the  case  of  other  neutron 
sources.  fbrJ.WUS  neutrons  v,=  42  keV  /tin.  and  lot 
6-MeY  monoenergetic  neutrons  q  =  18-7  keV  /on  in 
an  8-/tm  site,  which  is  similar  to  r,  Ibr  TRKiA 
neutrons  M.  (Jcrsienbrrg.  manuscript  in  prep¬ 
aration  .  whereas  i  |  =  I  or  15  7kcY/mi  were,  re¬ 
spectively.  assumed  in  the  model  calculations 
f Brenner  and  Hall  199(1  ;  thus,  depending  on  neu¬ 


tron  source,  the  model  input  variable  q  was  underesti¬ 
mated  by  up  to  a  factor  of  40.  so  it  follows  that  the 
magnitude  of  the  inverse  dose-rate  Ibr  these  sources 
was  overestimated  by  correspondingly  similar  factors. 
Applied  to  our  present  data,  this  means  that,  based 
on  the  Brenner  and  Hall  model,  no  enhancement  can 
be  expected  in  exponentially  growing  cultures  of 
C3H  10 1  I  cells  even  at  dose  rates  <  O  ')  cOy/min  in 
the  dose  range  examined.  To  further  illustrate  the 
impact  of  the  higher  q  on  this  model,  we  note  that  no 
dose-rate  offer l  ran  be  predicted  to  occur  in  neutron 
doses  and  dose-rates  relevant  to  radiation  protection. 
Specifically,  to  use  the  example  in  Brenner  and  Hall 
1 1990  .  no  enhancement  for  fission  neutrons  can  be 
expected  at  0-5  cGy  and  a  dose-rate  =  10"* cOy/min. 
In  fact,  we  calculated  that  in  order  to  obtain  the  two¬ 
fold  transformation  enhancement  for  fission  neutrons 
at  this  dose-rate,  the  total  dose  would  have  to  be 
50eGy.  Assuming  a  quality  factor  of  10.  this  corre¬ 
sponds  to  the  dose  which  is  100  times  the  annual 
effective  occupational  dose  permitted  by  current  I  S 
standards. 

In  summary,  it  seems  reasonable  to  conclude  from 
the  present  direct  comparisons  of  cell  survival,  trans¬ 
formation  and  age  distributions  of  control  and  irra¬ 
diated  cells,  as  well  as  by  examining  assumptions  in 
biophysical  model-  of  the  inverse  dose-rate  effect, 
that  there  is  no  real,  in  a  statistical  sense,  modifica¬ 
tion  of  these  end  points  by  extending  the  duration  of 
fission  neutron  exposure  or.  equivalently,  by  lower¬ 
ing  the  fission  neutron  dose-rate. 
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Appendix 

The  following  methods  are  independent  of  any 
model  assumptions  for  error  calculation  of  surviving 
fractions,  transformation  frequency  or  derived 
qualities. 


A1 .  Test  of  independence 

One  criterion  of  independence  is  that  the  mean 
DRMF  is  the  same  for  each  value  of  dose  or  dose- 
rate.  which  in  the  case  of  linear  regression,  requires 
that  the  slope  equals  0.  To  test  for  independence  of 
DRMF,  we  calculated: 

(b— 0)  •  SxyJn—  1 

t=  f  > 

•V* 

where  b  is  the  slope,  n  the  number  of  observations, 
and  Sx  and  Sy.x  have  the  usual  meaning  (Dixon  and 
Mascy  1957;.  and  compared  its  value  with  the 
/-distribution.  The  null  hypothesis:  b  —  0  can  be 
rejected  at  the  2-lcvcl  of  significance,  and  the  two 
variables  are  said  to  be  dependent,  if  t<tll2.x(n-2) 
or  if />/,  _  i  2 . , ! «  —  2i . 


A2.  Multiple-regression  analysis 

This  section  will  provide  explicit  formulas  and 
mathematical  procedures  used  in  the  multiple- 
regression  analysis  of  transformation  data  graphed 
in  Figure  3.  They  arc  also  suitable  for  a  similar 
analysis  of  any  results,  if  the  assumption  can  be  made 
that  all  the  curv  es  have  a  coinciding  intercept:  in  the 
present  case,  this  is  a  reasonable  expectation  since 
the  intercept,  by  definition,  is  the  spontaneous  trans¬ 
formation  frequency. 

We  consider  our  exposure-time  groups,  n/=0,  1, 
2.  3.  4  .  where  j  =  0  denotes  an  arbitrary  reference 
group,  chosen  here  to  be  the  high  dose-rate  group, 
and  separate  independent  variables  xk(k ~  1, 2.  3.  4;, 
such  that  xk—  1,  if  group  Ay  and  =0  otherwise.  The 
multiple-regression  equation  or 

‘full'  model:  T~  a  +  d-  b0  +  xx  •  d  ■  />,  4 - 1 -xA-  d-  bA 

has  six  independent  multiple-regression  coefficients 
to  be  found.  A  regression  coefficient  such  as  bx, 
denotes  the  regression  coefficient  of  T  on  variable  .v, 
that  one  would  expect  to  obtain  if  other  variables  in 
the  'full'  model  equation  x2,  x3  and  xA  had  been  held 
constant  experimentally.  Thus,  the  ‘full’  model  for 
group  I:  T=a  +  d ■  b0  + xx  ■  d ■  bx,  for,  group  2:  T 
=  a  +  d ■  b0  +  x2  -  d-  b2,  etc.  The  null-hypothesis  is. 


H0 :  bx  —  b2  —  b3  —  bA  —  0, 
and,  accordingly: 

H0  model:  T=a  +  d-  b0. 

It  should  be  noted  that  the  above  equation  doses 
not  imply  that  the  actual  relationship  between  trans¬ 
formation  frequency  and  the  dose  is  linear,  but 
rather  than  a  linear  approximation  explains  a  sig¬ 
nificant  amount  of  variability  in  T  over  the  range  of 
d  values,  which  is  true  in  the  case  of  our  data. 

The  procedure  compares  the  calculated  F  to  the 
Fa(/i1(/i2)  distribution  with  the  appropriate 
numbers  of  degrees  of  freedom;  the  symbols  have  the 
following  definitions:  =  [n  —  2)  —  (n  —  6),  p2  =  n  —  6 
are  degrees  of  freedom  in  the  numerator  and 
denominator,  respectively,  a  =  the  cumulative  prob¬ 
ability,  n  =  the  number  of  observations  in  all  groups, 
and 

[SS£(H0  model) — SSE  (‘full'  model)]//i, 

SSE  (‘full'  model) /p2 

where  S££=sum  of  squared  sampling  errors,  defined 
just  below. 

Let  Ijj  =  observed  transformation  response  of  ith 
group  atjth  dose,  and  dj—jth  dose. 

For  H0  model:  SSE=l[tij— (a  + dj- b0))2, 

where  summation  is  performed  over  groups, 
i(  =  0,  1, . ..,  4)  and  data  points  in  the  combined  set, 
j ( =  1 .  . . . ,  22  for  our  data  set:  Table  1 ). 

For  ‘full'  model:  SSE=  TT  T  -  TT  X  •  /?, 

where  the  matrix  X  is  defined  below,  T  denotes  a 
vector  constructed  from  /y  values,  and  /?  denotes  a 
vector  constructed  from  multiple-regression  coeffi¬ 
cients  found  by  solving  the  following  matrix 
equation: 

r~  — i  r .  —i 


lOA 

1  aj  0  0  0  0* 

l0.2 

1  d2  .0  0  0  0 

*1.1 

1  d6  10  0  0 

a 

^0 

b, 

bi 

bs 

Ua 

1  d20  0  0  0  1 

Ja. 

U.2 

1  d2i  0  0  0  1 

J*.3_ 

J  d22  0  0  0  1_ 

T  X  f{ 
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This  equation  has  the  solution: 

0=(XT  XrVX1  T; 

which  is  then  used  to  determine  SSF.  for  ‘full'  model 
and,  finally  F.  For  the  present  data:  b0=  -0-433. 
bi  =0-103, '*2  =  0-222,  i3= -0-363.  b4  =  0-354,  and 
i5=  -0-981.  The  sign  change  in  multiple-regression 
coefficients  indicates  that  both  up  and  down  adjust¬ 
ments  in  groups  1  (r  =  0-5h;,  2  (t=  I  h;.  3  it  =  3  hi. 
4  (t  —  4-5  h '  are  required,  consistent  with  small 
/’’-values  which,  in  turn,  corresponds  to  a  high 
probability  that  the  same  curve  holds  for  the  data. 
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